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Abstract 

This report proposes a new method to study vortex tubes in one-dimensional 

velocity data of experimental turbulence. Vortex tubes are detected as local 

maxima on the scale-space plot of wavelet transforms of the velocity data. 

Then it is possible to extract a typical velocity pattern. The result for data 

obtained in a wind tunnel is consistent with those of three-dimensional direct 

numerical simulations. 

PACS numbers: 47.27.Gs, 07.05.Kf 
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By using direct numerical simulations and bubble/ cavitation experiments, it has been 
established that turbulence contains vortex tubes Regions of intense vorticity are or- 

ganized into tubes. Their radii and lengths are, respectively, of the orders of the Kolmogorov 
length rj and the integral length L. Vortex tubes occupy a small fraction of the volume and 
are embedded in a background flow, which is random and of large scales. 

Especially when the Reynolds number Rex is high, effects of vortex tubes on the velocity 
field are of interest @] . The velocity signal at small scales is enhanced only in a fraction of the 
volume. This fine-structure intermittency is attributable to vortex tubes. With increasing 
the Reynolds number, the flatness factor of the velocity derivative becomes high. This means 
that turbulence becomes more intermittent. 

However, in direct numerical simulations, the Reynolds number is always low, Re\ ^ 200. 
For the higher Reynolds numbers, one inevitably resorts to standard experiments, where a 
measurement is made with a probe suspended in the flow, and merely a one-dimensional 
cut of the velocity field is obtained. To study vortex tubes in such one-dimensional velocity 
data, we propose a new method. 

The experiments often deal with the velocity components in the mean-flow direction 
alone (hereafter, the streamwise velocity u), but the transverse velocity v is more suited to 
detecting rotational flows such as those associated with vortex tubes We accordingly 

use this velocity component. 

Our method is based on orthonormal wavelets, i.e., self-similar functions localized both in 
scale and space. There are several known families of wavelets [^|10|- With a wavelet family, 
wavelet transforms are computed from the signal. Since a wavelet function has a zero mean, 
the transform corresponds to a variation in the signal at a given scale and a given position. 
Moreover, since a wavelet family constitutes a complete orthonormal basis, the transforms 
retain the same information as the original data. This advantage makes interpretation of 
the results reliable, especially in statistical analyses [0,|T^. 



Suppose that a function v{x) is sampled as a discrete signal v[n] = v{n6x), where 6x is 
the sampling interval (n = to 2^ — 1). The wavelet transformation is written as 



N-1 V~\ 

v\n\ = (t;) + ^ ^hkWi,kH. (1) 

j=0 A:=0 

Here (■) denotes the average. The integers j and k specify the scale of the wavelet function 
Wj,fc[n] and its position on the x-axis, respectively. The wavelet transform -Oj^fc is the inner 
product of v\n\ with the wavelet function Wj,fc[?T-] and represents a signal variation of the 
scale I = 2^~^~^5x around the position x = 2^~^k6x |TT|,p!2 . 



We systematically detect vortex tubes with different scales and strengths as local maxima 
on the scale-space plot of Vj^k/i'^j.k)- Here {Vj^k) is the second-order moment computed 
for each of the scales, corresponds to its mean energy, and reflects both tubes and the 
background flow. The circulation flow of a tube of the size i at the spatial position x is 
expected to enhance the value of Vj^k/ {'^j,k) ^^e scale i around the position x If 
the value of Vj^i^/ a certain scale-space position (j, k) is greater than those at adjacent 
positions, which are indicated by shaded areas in Fig. |I], we consider that a tube exists at 
that position. 

The present analysis is based on Haar's wavelets, each of which is a sharp pulse in space: 



-2'^ for n = 2^-^k to 2^-^k + 2^-^-^ - 1 

+2-^ for n = 2^-^k + 2^-^"^ to 2^-\k + 1) - 1 (2) 

elsewhere. 



Because of poor localization in scale, Haar's wavelets were not favored in previous studies. 
However, they have advantages. First, contrasting to other orthonormal wavelets that oscil- 
late many times in space, Haar's wavelets represent single oscillations. Since the transverse- 
velocity profile of a vortex tube has the same character (see below), Haar's wavelets work 
better in analyses of vortex tubes. With the other wavelets, we have analyzed our data. A 
tube is detected repeatedly at distances from the wavelet center. Second, a Haar's transform 
is analogous to a velocity increment 5v{t) = v{x + i) — v{x), the standard tool to analyze 
experimental turbulence [|12|. They have the same sign. Statistical results for Haar's trans- 



forms can be interpreted as those for velocity increments. 
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The most familiar model for a vortex tube is a Burgers vortex. This is an axisymmetric 
steady flow in a strain fleld. In cylindrical coordinates, they are written as 

2z/ 



Me oc — 
ar 



1 — exp 



2 ^ 

ar 



4z/ 



(a > 0), (3) 



and 



{Ur, Me, u^) = (^-^«^7 0, • (4) 

Here u is the kinematic viscosity. The above equation (|]) describes a rigid-body rotation 
for small radii, and a circulation decaying in radius for large radii. The velocity is maximal 

ro is regarded as the tube radius. 
Suppose that the vortex tube penetrates the (x, |/)-plane at the point (0, A). Here the 
X- and y-axes are, respectively, in the streamwise and transverse directions. If the direction 
of the tube axis is {6, (p) in spherical coordinates, the streamwise (w) and transverse (f ) 
components of the circulation flow Me are 

Acos6' xcos^ , . 
u = ue[r) and v = ue[r), (5j 

with 

= — sin^ 9 cos^ ip) + A^(l — sin^ 9 sin^ ip) + 2xA sin^ 9 sin ip cos v?. (6) 

Likewise, for the radial inflow Ur of the strain fleld (P, the streamwise and transverse 
components are 

x(l — sin^ 6^ cos^ 09) + A sin^ 6^ sino9 COS09 
u = Ur{r) [7) 

and 

a; sin^ 6' sin 09 COS09 + A(l — sin^ 6' sin^ oj) , 
V = Ur{r). (8) 

If a tube passes close to the probe (A < rg) and the tube is not heavily inclined (6* ~ 0), 
the transverse velocity is dominated by the circulation flow. This situation is important. If 



A ^ tq, the velocity signal of the tube is weak and diluted by the background flow. If ^ 0, 
the contribution of the circulation flow is small. The transverse velocity is then dominated 
by the radial inflow, which does not have a speciflc scale. To such scale-free motions, our 
method to detect tubes is insensitive. 

The velocity proflles of vortex tubes with A ^ tq and ^ ~ are nearly the same []7|,p!3 



However, as indicated by the flrst term of Eq. (P), the radius of a tube observed along the 
streamwise direction is different from its true radius ro as ^ ro/(l — sin^ 9 cos^ ^Y^"^- Thus 
self-similar functions such as Haar's wavelets capture efficiently those tubes. Although there 
might be contamination from tubes with A ^ ro or ^ 0, such a drawback is inherent in 
any of experimental works. 

Hereafter we analyze data obtained in a wind tunnel Its test section was 3 x 2 x 
18 m in size. Turbulence was produced by placing a grid across the entrance to the test 
section. The grid was made of two layers of uniformly spaced rods, the axes of which were 
perpendicular to each other. The separation of the axes of the adjacent rods was 40 cm. 
The cross section of the rod was 6x6 cm. We simultaneously measured the streamwise 
{U + u) and transverse (v) velocities with a hotwire anemometer. The probe was positioned 
on the tunnel axis at 6 m downstream of the grid. The signal was low-pass filtered at 8 kHz 
and sampled digitally at 16 kHz. The entire length of the signal was 2 x 10'' points. 

The mean streamwise velocity U is 8.70 m s~^. The root-mean-square fiuctuations (n^)^/^ 
and (t;^)^/^ are 0.446 and 0.427 m s~^, respectively. Since the turbulence level {u'^)^^'^/U 
is less than 10%, we rely on the frozen-eddy hypothesis of Taylor, d/dt = —Ud/dx, which 
yields the integral length L of 17.2 cm, the Taylor microscale A of 0.858 cm, the Kolmogorov 
length 1] of 0.0270 cm, and the microscale Reynolds number Rex of 260. 

The velocity signal is divided into 2400 segments of 2^^ points. To each segment, we apply 
the wavelet transformation (A^ = 13). Then statistics are computed over the segments. We 
present results only for the scale i = 81] {j = 10), i.e., the smallest scale to which our method 
to detect vortex tubes is applicable. Results at the other scales are similar. 

Fig. 1^ shows probability density functions (PDFs) of Haar's transforms (solid lines) 



and velocity increments (dotted lines). Since the transverse- velocity PDFs are symmetric, 
those of the absolute values are shown. The agreement between the wavelet transforms and 
velocity increments is excellent. 

By averaging signals centered at the position where -Ol;./ is locally maximal, we 
extract typical patterns of vortex tubes in the streamwise (u) and transverse (v) velocities. 
The detection rate of the local maxima per the integral length L is 3.37 at £ = 877. Since 
a wavelet function is spatially extended, we determine in each case the center position so 
that the absolute value of the velocity increment is maximal. When the increment is 

negative, we invert the sign of the v signal. The result is shown in Fig. |^ (solid lines). The 
u pattern is shown separately for 6u > and 6u < (designated, respectively, as and 
u'). We also show velocity profiles (^) and (|^) of a Burgers vortex. It is assumed that the 
tube center passes through the probe position (A = 0) and the tube axis is perpendicular 
to the streamwise and transverse directions (6 = 0). The tube radius r^. is determined so as 
to reproduce the v pattern. 

The V pattern of grid turbulence is close to the profile of a Burgers vortex. Since the v 
pattern is somewhat extended, there might be additional contributions from vortex tubes 
with A ^ ro or ^ and vortex sheets. Previously, velocity patterns of tubes were studied 
by averaging for large values of wavelet transforms or velocity increments The pattern 

shapes are close to those in Fig. ^ but are biased toward strong tubes. Our results based 
on local maxima are more representative of vortex tubes in turbulence. 

The patterns of grid turbulence appear to be dominated by the circulation flow Me 



of a vortex tube [jT3[. There is no significant evidence for the presence of the radial inflow 
Ur- A vortex tube is not necessarily identical to a Burgers vortex. The same conclusion was 
obtained from direct numerical simulations 

The tube radius = 101] observed at i = 8rj serves as an upper limit for the intrinsic 
tube radius ro, which should be several of the Kolmogorov length rj. With an increase 
of the scale i, there is increasing importance of inclined vortex tubes. There also exist 



tubes with intrinsically large radii |]14|. On the other hand, the amplitude of the v pattern. 



corresponding to the circulation velocity of a typical tube, is of the order of (f^)^/^. The 
radius of a tube and its circulation velocity estimated here are consistent with those obtained 
from direct numerical simulations |jl|J^. 

The statistical fluctuation of the transverse velocity v is shown in Fig. ^. Vortex tubes 
are embedded in a large-scale background flow. Hence, from the velocity signals around the 
local maxima of v'jj^/ (v^k), the expansion formula (|1|) had been used to remove motions with 
scales i > 81] {j < 10). The resultant fluctuation, which mainly reflects differences of tube 
parameters, is comparable to that in a direct numerical simulation 0. 

The PDFs of wavelet transforms and velocity increments for the transverse velocity v in 
Fig. |^(b) have long tails 0,||,|l2|. It is possible to study this intermittency phenomenon at 
each scale by comparing a PDF of wavelet transforms at local maxima of Vjk/{''^jk) "with 
a PDF of the whole transforms as in Fig. |] (solid lines). The PDF of the local-maximum 
transforms has a peak at l-Oj^fcl ~ {Vjk)^^'^ ^^^^ ^ long tail toward large magnitudes. If the 
condition for the local maximum is relaxed so as to ignore transforms at the smaller scale 
j + 1, the result reproduces perfectly the tail of the whole-transform PDF (dotted line). 
Thus vortex tubes account for the long tail of the PDF. The situation is the same in veloc- 
ity increments, since Haar's transforms are statistically equivalent to velocity increments. 
However, only 10% of the scale energy is shared by the above subset of wavelet transforms 
that reproduce the tail of the whole-transform PDF. The background flow is energetically 
predominant |]T^ . 



Finally, we give additional comments. First, our method to study vortex tubes is crude. 
They have been studied with three-dimensional direct numerical simulations in more reli- 
able ways HIJ^. Nevertheless, our method is useful at high Reynolds numbers, where only 
one-dimensional experimental data are available. Second, although our method is based on 
wavelet transforms for the transverse velocity alone, those for the streamwise velocity are 
also useful. For example, they could constrain the local strain field (see Fig. 0(a)). Third, 
our method is intermediate between those of our previous works. We proposed statistical 
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measures based on orthonormal wavelets |T2[. Although these measures are not relevant 
directly to tubes, they are robust with respect to the choice of the wavelets and thus char- 
acterize rigorously scale-space structures of a velocity signal. We also proposed a method to 
detect tubes by using the transverse-velocity profile of a model tube as a nonorthonormal 
wavelet Although a nonorthonormal wavelet transformation is redundant, this method 
provides high-resolution estimates of the size and spatial position of a tube. We hope that 
the application of these methods would improve knowledge of small-scale structures of tur- 
bulence. 
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FIGURES 

FIG. 1. Schematic representation of wavelet transforms on a scale-space plot. If the value of 

Vj j^/{v'j f^) at the position (j, k) exceeds tliosc at the adjacent positions (shaded areas), we consider 
that v^k/iv^k) is locally maximal at (j, fc). 

FIG. 2. PDFs of wavelet transforms (solid lines) and velocity increments (dotted lines) at 
£ = 8r] for the streamwise (a) and transverse (b) velocities. The abscissa is normalized by the 
root-mean-square value: = 0-124, {Su'^y/^ = 0.0950, = 0.169, and {6v^)'^^'^ = 

0.132 m s~^. We also show Gaussian distributions with zero means and unity standard deviations. 

FIG. 3. Conditional averages of the streamwise (u) and transverse {v) velocities for local 
maxima of v^jk/{v^^ at £ = Sr] (solid lines). The abscissa is the spatial position x normalized by 
the Kolmogorov length 77. We show the streamwise velocity separately for 5u > {u~^) and 5u <0 
{u~). The dotted lines represent profiles of a Burgers vortex. 

FIG. 4. Conditional average of the transverse velocity v for local maxima of v^k/iv^k) at 
£ = 877. The abscissa is the spatial position x normalized by the Kolmogorov length 77. The error 
bars denote the statistical fluctuation (it la), which has been computed after we had removed 
motions with scales £ > Srj. 

FIG. 5. PDFs of wavelet transforms of the transverse velocity \vj,k\ at ^ = 877. The abscissa 
is normalized by the root mean square of the whole transforms at the scale. We show 

PDFs for local maxima of Vj k/i'^jk) whole transforms (solid lines), normalized by the 

number of the whole transforms. The fraction of the scale energy shared by the local-maximum 
transforms is 5.24%. We also show a PDF for the local maxima determined without the smaller 
scale j + I (dotted line). Their energy fraction is 9.56%. 
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